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Recently [J. Haro and E. Elizalde, Phys. Rev. Lett. 97, 130401 (2006)], a Hamiltonian formulation has been 
introduced in order to address some longstanding severe problems associated with the physical description 
of the dynamical Casimir effect at all times while the mirrors are moving. Here we present the complete 
calculation providing precise details, in particular, of the regularization procedure, which is decisive for the 
correct derivation of physically meaningful quantities. A basic difference when comparing with the results 
previously obtained by other authors is the fact that the motion force derived in our approach contains a reactive 
term — proportional to the mirrors' acceleration. This is of the essence in order to obtain particles with a positive 
energy all the time during the oscillation of the mirrors — while always satisfying the energy conservation law. 
A careful analysis of the interrelations among the different results previously obtained in the literature is then 
carried out. For simplicity, the specific case of a neutral scalar field in one dimension, with one or two partially 
transmitting mirrors (a fundamental proviso for the regularization issue) is considered in more detail, but our 
general method is shown to be generalizable, without essential problems (Sect. 2 of this paper), to fields of any 
kind in two and higher dimensions. 

PACS numbers: 42.50.Lc, 03.70.+k, ll.lO.Ef 



I. INTRODUCTION 

Moving mirrors modify the structure of the quantum vacuum, what manifests in the creation and annihilation of par- 
ticles. Once the mirrors return to rest, a number of the produced particles will generically still remain, which can be 
interpreted as radiated particles. This flux has been calculated in the past in several situations by using different methods, 
as averaging over fast oscillations 1110], by multiple scale analysis |3], with the rotating wave approximation |4], with 
numerical techniques jsj], and others la]. Here we will be interested in the production of the particles and in their possible 
energy values, all the time while the mirrors are in movement. This is in no way a simple issue and a number of problems 
have recuiTently appeared in the literature when trying to deal with it. To start with, it is far from clear which is the 
appropriate regularization to be used. Different authors tend to use different regularizations, forgetting sometimes about 
the need to carry out a proper (physical) renormalization procedure, in order to obtain actually meaningful quantities. 
Thus, it turns out that ordinarily, in the case of a single, perfectly reflecting mirror, the number of produced particles as 
well as their energies diverge all the time while the mirrors move. Several prescriptions have been used in order to obtain 
a well-defined energy, however, for some trajectories this finite energy is not a positive quantity and cannot be identified 
with the energy of the produced particles (see e.g. |01)- 

Our approach relies on two very basic ingredients First, proper use of a Hamiltonian method and, second, the 
introduction of partially transmitting miiTors, which become transparent to very high frequencies. We will here prove in 
this way, both that the number of created particles remains finite and also that their energies are always positive, for the 
whole trajectories corresponding to the mirrors' displacement. We will also calculate from first principles the radiation- 
reaction force that acts on the mirrors owing to the emission and absorption of the particles, and which is related with 
the field's energy through the ordinary energy conservation law. This carries along, as a consequence, that the en ergy o f 
the field at any time t equals, with the o ppo site sign, the work performed by the reaction force up to this time t 1191 [lOll . 
Such force is usually split into two parts ill ll[T2ll : a dissipative force whose work equals minus the energy of the particles 
that remain |9], and a reactive force, which vanishes when the mirrors return to rest. We will also prove below that the 
radiation-reaction force calculated from the Hamiltonian approach for partially transmitting mirrors satisfies, at all time 
during the mirrors' oscillation, the energy conservation law and can naturally account for the creation of positive energy 
particles. Also, the dissipative part obtained within our procedure agrees with the one calculated by other methods, 
as using the Heisenberg picture or other effective Hamiltonians. Note, however, that those methods have traditionally 
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encountered problems with the reactive part, which in general yields a non-positive energy that cannot be considered as 
that of the particles created at any specific t. 

The organization of the paper is as follows. In Sect. 2 we introduce first the canonical formulation underlying the 
whole procedure. In particular, we give the explicit expressions for the Hamiltonian and the corresponding energy. We 
do this by considering the Hamiltonian method for a neutral Klein-Gordon field in a cavity, generically in (ri + 1)— 
[although we will mainly work in (3 + 1)—] dimensional space-time, with boundaries moving at a certain speed v ^ c. 
In Sect. 3 we deal with the case of a single partially transmitting mirror. We formulate the quantum theory based in the 
Hamiltonian approach and, successively, the quantum theory in the Heisenberg picture. We finish that section with a 
detailed comparison with early known results, e.g. those obtained with the method of Jaekel and Reynaud, and with the 
method of Barton and Calogeracos. In Sect. 4 we study the more difficult case of two partially transmitting mirrors. In this 
situation, the part of the Hamiltonian that describes the interaction between the field and the mirrors depends on e what, 
in general, renders it quite difficult to describes this part. For that reason, the reactive part of the motion force can seldom 
be calculated. In any case, we prove here, in particular, that our dissipative part of the motion force exactly coincides 
with the dissipative force obtained by Jaekel and Reynaud ifisll . Moreover, following the Hamiltonian approach we show 
that the problem of the negative energy that appears in the Davis-Fulling model can be resolved if partially transmitting 
mirrors are considered, what is in our view a very physical approach to the renormalization issue. The last section of the 
paper is devoted to a final discussion and conclusions. 



II. CANONICAL FORMULATION OF THE PROBLEM 



We consider a neutral massless scalar field in a cavity, Q,t, and assume that the cavity boundary is at rest for all times 
t < and returns to its initial position at time t = T, to remain there for a while. Its velocity will be given in terms of c, so 
that we will work with the dimensionless quantity e — v/c ^ 1. In a practical situation, as the one featured in Ref. Iil4i1 . 
this turns out to be of order 10^* (more about that later). 

The Lagrangian density of the field is 

C{t, x) - i [{dt<f>f - |Vx<^|'] , Vx e r!t c M", Vi e M, (1) 
If we use the canonical conjugated momentum 

dC 

the energy density of the field is given by the expression 

£{t, x) = ^dt^ ~ C{t, x) = i {e + |Vx0p) , (3) 

and the energy itself is 

E{t-e)= [ d"x£{t,x). (4) 



A. Hamiltonian and energy 

It is a well-known fact that the energy density does not generically coincide with the Hamiltonian de nsity ifTsll - lfTTll . 
Here, to obtain the Hamiltonian density of the field we follow the method discussed in Refs. U8i] and lfl9ll . First, we 
transform the moving boundary into a fixed one by performing a (non-conformal) change of coordinates 

n : {i, y) ^ (t(t, y), x(i" y)) - [t, R(t" y)), (5) 

that transform the domain fit into a domain ft which is independent of time. Making use of the coordinates (t, y), the 
action of the system behaves as 



f f d^ydiC{i,y), 
Jm Jn 



(6) 
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wi\h C{t,y) = J£ (7?. (t, y)), where we have introduced the Jacobian J of the coordinate change, defined by (i"x = Jd^y. 
Let us now consider the function (j) given as 0(t, y) = \/j (f){'R{t, y)). Then, the canonical conjugated momentum is 

= ^/J9t0(7^(^,y)), (7) 
and, from here, the Hamiltonian density is obtained as 

nil y) = ^5,-0 - C{1 y) - i (e' + J| Vx0|') + ^(5*^ - ^^0). (8) 
In the coordinates (i, x), the Hamihonian density is given by 

H(i,x) = 7^(7^-l(<,x))^ - i7^(7^-l(^,x)). (9) 
Now, from expressions dlJ and O we have that ^(f, y) = \/j£_{TZ{i, y)), and a straightforward calculation yields 

7i(t,x) = £(i,x) +e(t,x) < a,R(7^-l(^,x)), Vx<^(t,x) > +i e(i,x)0(i,x)9,(ln J)|^_i(, ,^) . (10) 

B. A simple and explicit example 

As a simple example, in the case of a single mirror following a prescribed trajectory (t, eg{t)) in a 1 + 1 space-time, 
we can take R{t, y) ^ y + £g{t), and thus we explicitly get 

nit,x) ^ £it,x) + eg{t)^it,x)d^(l){t,x). (11) 

III. SINGLE, PARTIALLY TRANSMITTING MIRROR 

In this section we consider a single mirror in 1 + 1 space-time following a prescribed trajectory {t, eg{t)). When the 
mirror is at rest, scattering of the field is described with the 5*— matrix 

^^''> - [ r(^)e2-i sico) ) ' ^^^^ 

where x — L is the position of the mirror. The matrix S is supposed to be real in the temporal domain, as well as causal, 
unitary, and transparent to high frequencies fidl. More specifically, these conditions appear naturally as a consequence of 
the following considerations. 

1 . Since the field is neutral, it should be: 

S{-uj) = S*{oj). (13) 
In fact, the quantum field, in the Schrodinger picture, can be decomposed as 



where 

g^,R{x) = -L= {s{Lu)e-'^-9{L - x) + (e"^'^^ + r(c^)e-2-^e'--) 0{x - L)} , (14) 
g^,L{x) = {(e^"" + r(w)e2''"^e-*"") e{L - x) + s(u;)e*""0(i - x)] , (15) 
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are the right and the left incident modes, respectively 111211 . As is the usual procedure in Quantum Field Theory, 
when uj < one performs the change (iuij a^^ j, and thus, the field behaves as 



-B T.Jo ^ ^ 



(p[X) = 

j=R,L 

Now, since the field is neutral, it follows that 'g^^,j{x) — j{x), and finally, we conclude that 

S{-oj) = S*{uj). 

This proves the statement. 

2. As a consequence of the commutation rule [(j){t, x), (p{t, y)] = 0, it follows that 

S{uj)S\uo) = Id. (16) 
This is straightforward and needs no further comment. 

3. And, as a consequence of the commutation rule [^(f, x), y)] = ~iS{x — y), we obtain the following causality 
condition 

[ dujr{w)e-"^' - ' -J- ■ 

JR 

in a distribution sense, i.e.. 



duj r{w)e-"^' / dij s{w)e'"^' = 0, Wt<0, 



lim / duj r{w)p^{uj)e = Um / du! s{w)py{uj)e 



lim / dujr(w)p^(uj)e-''^* = lim / du s(w)p^(Lo)e-"^* = 0, Vt < 0, 

7- 

where is a frequency cut-off. This condition is satisfied when 

S{ui) is analytic for Im(a;) > 0, and s and r are meromorphic functions. (17) 

4. A physical mirror is always transparent to high-enough incident frequencies, thus it must necessarily hold 

^ Id, when |w| -> oo. (18) 

A. Quantum theory based on the Hamiltonian approach 

In order to obtain the quantum theory, we will work in the coordinates defined in the example above. In those coordi- 
nates the mirror is situated at the point y = and the right and left incident modes are given by Eqs. ( fT4] l and (flST l. with 
L = Q. Then, in the coordinates {t, x), the instantaneous set of right and the left incident eigenfunctions, which generalize 
the set used in the case of a perfectly reflecting mirror, is 



{t,x;e) = g^,,{x^eg{t)) j^R,L. (19) 



Note that, in general, we do not know explicitly the part of the Hamiltonian that describes the interaction between the field 
and the mirror. As a consequence, in order to obtain the quantum theory, the energy of the field E{t) = dx£{t, x), 
which in the presence of a single mirror does not depend on e, should be viewed as part to the free Hamiltonian of the 
system. 

As is usual, working in the interaction picture, the field is expanded as follows: 

/■OC 

^i{t,x;e)^ / dujaw,je^'"^^guj.j{t,x;e) + he, 



j=R,L 







/■oo 

ii{t,x;e) ^ -i ^ / dujujaw,je^''^'''gujj{t,x;e) + he, (20) 
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where he means, in each case, the Hermitian conjugate of the preceding expression. The quantum equation, in this picture, 
is given by 



2 



dx£,i{t, x; e)dx4>i(t, x; e) + he 
dx£,i{t, x; 0)dx4>i{t, x; 0) + he 



m 



(21) 



Let now T* be the quantum evolution operator of the Schrodinger equation (pT|) . and let |0) be the initial quantum state. 
Then, the average number of produced particles at time t is 



j=R,L • 

and the dynamical energy at time t, that is the energy of the created particles at time t, is obtained as 



(22) 



(23) 



j=R,L 



We should note that, during the movement of the mirror, the particles have been called sometimes quasi-particles, owing 
in part to the difficulties encountered in the past when trying to give them a physical sense (see ll2lll ). 
A simple but rather cumbersome calculation yields the following results 



(w + ^O' 



(24) 



m)) 



e 

4^ 



oc />oo 



dtoduj'ujuj' 

LU + Lu' 



(25) 



where 9t is the Heaviside step function at point t, e.g., 0((t) = 6{t — t), and f{uj) = j^dr f {T)e'^'^'^ is the Fourier 
transformed of the function /. These two quantities are in general convergent. However for the Davis-Fulling model |7] 
— that is, in the case of a single perfectly reflecting mirror — such quantities are divergent when the muTor moves or when 
the displacement exhibits some type of discontinuity lfT7ll22ll . 

In this situation, in order to obtain a finite energy, several authors have used different regularization techniques ll9ll- lfTlll . 
For example using a frequency cut-off e~'^'' with < 7 <C 1, the regularized energy is (see Ref. |il9J) 



{E{t-^)) = - 

OTT 



7r7 



mm 



f{T)dT 



Thus, imposing the kinetic energy of the moving boundary to be 

1 



exp 



(26) 



(27) 



where Alexp is the experimental mass of the mirror, those authors conclude that the renormaUzed dynamical energy, 
namely Enit), is (see Refs. 101, jgll-lfTlll) 

{En{t)) EE |1 \~mm + f f{r)dT] + 0{e^). (28) 



e 

Gvr ^ 

However, when t < 5, with < 5 <C 1 such renormalized energy is negative. This shows that, when the mirror moves, 
the renormalized energy cannot be considered as the energy of the produced particles at time t (see also the paragraph 
immediately after Eq. (4.5) in Ref. |7]). 

From our viewpoint, such meaningless result is just due to the fact that a perfect reflecting mirror is used in the deriva- 
tion, what is not physically feasible at any price. Physical mirrors will always obey a transparency condition of the kind 
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here proposed ( flST l. and then it comes out for free that the average number of produced particles and the dynamical energy 
turn out to be well defined and are both positive quantities. 

We have also calculated the radiation-reaction force. For a single mirror this force is the difference between the energy 
density of the evolved vacuum state on both sides of the mirror, namely, 

{pHait)) ^ lim ((0| (r*)t£(t,65(t) - |<5|)r*|0) - (0| {T')^ E{t,eg{t) + \5\)T'\Q)) , (29) 

where the subindex Ha means that the radiation-reaction force has been calculated in the Hamiltonian approach. We 
obtain 

/A / e f°° dujduj'ujuj' „ r , i , , 

Stt^ Jq Jq uj + uj' I J 

x(|rM + r*(c.')|' + - s*{J)\^) + 0{^). (30) 

An important remark is here in order Note that, as a consequence of the energy conservation law, the dynamical energy at 
time t is equal to minus the work performed by the radiation-reaction force up to time t (see Refs. |9] and 1 10]). This law 
is naturally satisfied if we use the Hamiltonian approach. It is then clear that dZSl l and dSOl l are related through the formula 

{E{t)) - -e / {FHa{T))g{T)dT. (31) 

Jo 



B. Quantum theory in the Heisenberg picture 

Following the method of ll23ll . we have calculated the "in" modes when the mirror describes the prescribed trajectory 
(i, eg{t)). Using light-like coordinates, u = t + x and v = t — x, the "in" modes can be written as 

«; 0) = { [siLo)e-''^^ - AUv; 0)] eieg{t) - x) 

+ [e-'^^ + r(cc;)e-*"" - B^{u; 0)] Oix ~ eg{t))) + ©(e^), (32) 
0) = { [e-'"" + r{Lo)e'^-^ + B^{v- 0)] d{eg{t) - x) 

+ [s{u;)e-''^^ + A^{u; 0)] e{x - eg{t))} + 0{e^), (33) 

where 



BUy;i) = ^ / dLo'e^-'ygi-u -u')[r*iu') +riLo)]e-''\'^'\. 
The average number of produced particles after the mirror returns to rest is 



,, „_D T -JO Jo 



|2 



i.j=R,L 



(34) 



where {F,G) = i J^dx{F*dtG — GdtF*). To obtain an explicit result, we calculate the Bogoliubov coefficients 
{(/)°"|, f/*™*, } in the null future infinity I+, because the "out"-going modes are a very simple expression in I+. The 
final result is (see Ref. lO): 



^2 poc foo 



Jo 



M{t>T)^—f I dLjduj'LJu;'\g{Lj + Lj')r{\r{u;)+r*{cj')\^ + \siu;)~s*{Lu')\^)+0{e^). (35) 
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From this expression it is not difficuh to calculate the number of particles at time t. We only need to consider the function 

' ~ 1 git), when s > t, 

because ^((w + lo') = ^^^i 9^t{'-^ + w'). Then, inserting this expression into ((35l) we obtain formula (|24)) . 

The radiation-reaction force calculated in the Heisenberg picture, namely {Fh [t)), is the difference between the energy 
density of the "in" vacuum state on the left and on the right side of the mirror. A simple calculation shows that the energy 
density of the "in" vacuum i24ll 

/•oo 

x))^ / (5«C.("' ^; ^)du<t^::M^ 0) + ^; o)9„0™*(^, v- o)) , 

on both sides of the mirror, is 

{£{t,x)) = / dajw ± — ^ / dujduj'ujuj' g{uj + oj')x{^^) 
Jo on Jr2 

X (1 + r{u;)r{u;') - s(cj)s(cj'))e"*''^+"'^"6l(±(e(7(t) - a;)) + ©(e^), (36) 

where x(^) = ^(^) ~ 6'(— is the sign function. Note that the term of order e is ill-defined, because the function 
ujLd'giu) + Lo'){\ + r{uj)r{uj') — s{ijj)s{lo')) is not integrable in and, to obtain a well-defined quantity, appropriate 
regularization is needed. 

If we define the regularized energy by 

{£{t, X- I '^^^''^ j)du^t:*{u, V- 7) + v; i)d,4>Z^{u, v; 7)) , (37) 

then the regularized motion force, in the Heisenberg picture results 



{FH{t; 7)) = — / dujduj'uJLu'g{Lu + uj'){x{uj) + x(^')) 



8^ 

X (1 + r{u;)r{Lu') - s(tj)s(tj'))e"^'l"l+l'^'l^e-*("+"''* + ©(e^). (38) 



This integral is convergent and already cut-off independent. Thus, a natural definition of the renormalized radiation- 
reaction force is 

{FH,ren{t)) = / dtoduj' WLo' g{uj + uj'){x{lo) + x{^')) 

X (1 + r(tj)r(w') - s{uj)s{u'))e-'^'^+'^"^' + 0{e^). (39) 



Two more remarks are in order. In general, expression (p9|) disagrees with the radiation-reaction force pO)) obtained 
using the Hamiltonian approach. A detailed discussion of this point will be given in the next section. 

It is also important to stress that other definitions of the regularized energy density give rise to different motion forces. 
It is, of course, also possible to obtain the radiation-reaction force (l30l l within a specific, tailored regularization. To 
prove such statement let us consider, for a moment, the usual case in the literature of a perfectly reflecting mirror. In the 
Heisenberg picture the annihilation operators on the left side of the mirror are 



dio' g9t{uj + cj')a^, 

LOU 



\/ bJLj' 1 
TT UJ + ljj' 



daj'g9t{uj - u!')a^'^L V 



where V denotes Cauchy's principal value. Then, on the left side of the mirror, the field — in the Heisenberg picture- 
can be written as follows 

4>Hit,x]e) = / duj(^d^^L{t)gi^,L{t,x;e) +dlj^{t)g'^j^{t,x]e)j 



8 



and, after some algebra, we get 
where 

fuj,L{t,x;j) = — [sin(a;a;) - eg{t) cos(a;a;)] 



Using the formula 



we easily find that 



-'-^ I dij'dTe'^'*^t{-uJ-uj')V ( -^—] sin(w'a;)e-^'" 
aui P e ^ ' — me ^ x(u — t), 



cj^:\{t,x-Q) = ^sm{ujx)e-''^' - 2ieJ —e-'^^^-'^^git - x) ^ ^^,L{t,x;0), 

and making the changes e ^ — e and x —x, we obtain the expression for the right modes. That is, we have obtained 
an equivalent expression as for the "in" modes. 

Then, defining the regularized energy from this new expression of the "in" modes 

1 f°° 

{£{t,x;-f)) = - ^ / duje~''''^{dt(pu,j{t,x;j)dt(plj{t,x;j)+d^ip^j{t,x;-f)d^(plj{t,x]j)), 
we obtain the following regularized motion force 



I" Jo Jo ^ + ^' 

which agrees with ((30l) . for the case of a perfectly reflecting mirror. 

The fact that different (sometimes a priori quite reasonable) regularization procedures may lead to very different finite 
results is well known |25], as also the ensuing consequence that there is generically no physics associated with an arbitrary 
regularization prescription, which calls for a subsequent renormalization procedure to establish contact with the physical 
world. This is more so when one deals with plain mathematical and physically unrealistic boundary conditions, as we 
have learnt from a number of situations involving the ordinary Casimir effect too |26]. However, this essential point 
seems to have been put aside, at least to some extent, when dealing with the problem at hand, maybe due to the intrinsic 
mathematical difficulty of this issue here. This was the motivation for the last explicit calculation above, that we consider 
a sufficiently clarifying exercise which exhibits what is going on here. In the following we will proceed with a strict 
comparison of our results with those of other authors that have previously appeared in the Uterature on the subject. 

C. Comparison with Isnown results 

(i) The method ofJaekel and Reynaud. 

To study the radiation-reaction force, these authors ll27ll consider the following effective Hamiltonian 

Hj^R = -eg{t)F{t), (40) 



where {t, eg{t)) is the trajectory of the miiTor and F{t) = lim^^^o \£{t, ~ £{t, +|<5|)j is the force operator at the 
point a; = 0. A simple calculation yields 

/•oo poo 

{FjMt)) = I dcudiu'cuuj' Im [e-^(-+-')*5?t(c^ + uj') 

x(|rH + r*{uj')\' + \siu;) - s*{l,')\') + 0{e'). (41) 
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Integrating by parts, we obtain 

Jo Jo UJ + UJ' 

what shows that expression (|4T]) is divergent while the mirror moves. To obtain a regularized quantity we write (|4T]) as 
follows 

{Fj^R{t)) = ^ f drgir) f dojdoj'oju:' (x{uj) + x{^')) [(1 + r{uj)r{oj') ~ s{u:)s{u:')) 

+ (1 + r*{Lo)r*{uj') ~ s*{lo)s*{lo'))] + O(e^), 

and we define the regularized motion force by 

{Fj,R{t;-f)) = / dTg(r) / dwdw'ww'(x(w) +x(^')) [(1 + r(w)r(w') - s{iJ)s{Jy)p^{ijj,J) 

J-oo JR2 

+ (1 + v*{loY{J) - e-^^^+'^'^f*-^) + 0(e2), (42) 

where the cut-off p.^(yi-, lo') is a meromorphic function, analytic for Im(a;) > and Im(a;') > 0. 

Now, applying the causality of the S'-matrix (Eq. (fTTIi). and making 7 — > 0, an easy calculation leads us to the expression 
(|39] l. For this reason, defining the renormalized radiation-reaction force through the formula (|39[). i.e., 

{F.J,R,ren{t)) = {FH,ren{t)), (43) 

one does conclude that the method of Jaekel and Reynaud is equivalent to the quantum theory in the Heisenberg picture. 
Note also that 

dt{FHa{t))m^^ / dt{Fj 

,R,ren 

Jr 

This identity proves that the dissipative parts of (Fffa(i)) and {Fj,R,ren{t)) always agree. 

On the other hand, in several situations, the reactive part disagrees actually. For example, if r{w) = — ^!|° ^ , and 
s(w) — f- , with a > 0, there holds the relation 

{FHa{t)) = + {Fj,R,ren{t)), (44) 

where 

{Fj,R,ren{t)) = ^ dz dr (j^ - g-^^^*--) 5 (t). (45) 

That is, both motion forces differ in a reactive term. Note also that, during the oscillation of the mirror, the work done 
by the motion force (Fj,fl,ren(0) a negative quantity. Consequently, from the previous remark it follows that 

the dynamical energy is not positive, and therefore a seemingly meaningless result is obtained since, in our opinion, the 
dynamical energy is to be interpreted as the energy carried out by the produced particle. To avoid such difficulty the 
reactive term ~§^ 'g{t) should not be arbitrarily suppressed but, on the contrary, has to be duely taken into account. This 
saves the day and endows the whole picture with physical sense, as explained in the previous section. 

(ii) The method of Barton and Calogeracos. 

In Ref. |fl5] (see also |28] and |29]), these authors study the particular case r{'w) = — ^^^^ , and s{w) — ^ with 
a > 0. In such situation, the interaction between the field and the mirror can be described by the Lagrangian density 

i [{dtcl))^ - (9,0)2] - a(l>^{x ~ eg{t)). (46) 



Following the method discussed in Sect. II, we have obtained the quantum Hamiltonian 



H(t)^ / dxE{t,x)+a(j)\t,tg{t)) + 



dx ^ (i, x)dx4'{t, x) + he 



(47) 
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Now, inserting (|20] | into the integral dx£{t, x) + acfp' {t, eg{t)), we get 




We thus conclude that the quantum equation in the interaction picture is given by expression ( 1211 1. And, consequently, for 
these reflection and transmission coefficients, the authors would obtain the same formulae jTM . ( |25] l and dSOl ). 

It should be noted, however, that two relevant differences exist between their results and the ones we have derived here 
previously. 

1. In order to obtain the quantum equation, B-C make a unitary transformation which does not seem to be easily 
generalizable to the case of two moving mirrors. In our case this can be done without the least problem, as we shall 
show below. 

2. In the above mentioned paper, Ref. fl^l, the authors use the same technique for mass renormalization that had been 
employed by Barton and Eberlein in Refs. 1 10] and 1 1 1] for the case of a completely reflecting mirror, in order to 
eliminate the reactive part of the motion force. However, within such renormalization, the energy of the field is not 
a positive quantity for all time t and, consequently, the concept of particle is again ill-defined during the oscillation 
of the mirror 



IV. TWO PARTIALLY TRANSMITTING MIRRORS 



In this section we consider the situation where we have two moving mirrors which follow prescribed trajectories 
{t, Lj{t; e)), where Lj{t;e) = Lj + egj{t), with j ~ 1,2, and we assume that Li{t;e) < L2{t;e), Vt E M. In this 
case it is impossible, in practice, to work within the Heisenberg picture, because it is indeed very complicated to ob- 
tain the "in" and "out" mode functions in the presence of the two moving mirrors. Alternatively, with the purpose to 
calculate the dissipative part of the motion force, the number of radiated particles, and their energy, one can use the 
approach due to Jaekel and Reynaud and based on the effective Hamiltonian Hj^r = 2 ^9ji^)-^ji^)' where 

Fj(t) = lima^o (^£{t, Lj — \5\) — Lj + \S\)^ is the force operator at the point x — Lj (see Refs. |27] and fl3\). 

However this method is not useful in order to calculate the reactive part of the motion force or the dynamical energy while 
the mirrors are in movement. 

To obtain those last quantities we are led to use, once more, the Hamiltonian approach. In this case, if we consider the 
change 



R{ly) 



the Hamiltonian density of the field is 



1 



L2 — L] 



[L2 {t; e)iy-L,) + Li{t; e) {L2 - y)] , 



j~i,2 L2{t;e)- Li{t;e) 



d^(j){t,x){x - Lj{t; e)) + -<j){t,x) 



(49) 



where L^i^^ {t; e) = Lpj (t; e). In the coordinates {t, y) the set of right and left incident modes can be obtained from 
Eqs. (8) and (9) of Ref. UJ]. We find 



1 I Sl{uj)s2{uo) 



S2{i^) 

d{uj) 



e-^^y0{L, - y) 



-twy 



d{uj) J 

e{y-L2)), 



d{Lu) 



(50) 
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d{Lu) d{uj) 

Si{uj)s2{uj) 

d{to) 



2iwL2 \ ~iwy 



V d{uj) 

9{y - Li)9{L2 - y) 

e'^y0{y~L2)^ 



0{Li-y) 



(51) 



where d{uj) = 1 — ri{uj)r2{uj)e^'^'^^-'^^ Then, the instantaneous set of right and left incident eigenfunctions in the 
coordinates (i, x) is 



L2 — Li 



r5<^,fl(2/(^, 2;)), 9oj,L{t,x;e) 



L2 — Li 



_ L2{t-e)- Li{t-e) 
The fields can by expanded as follows 

/■OO 

4){t,x)= ^ / dujdu,jgujjit,x;e) + he, 

^ / dLJUjd^^jg^j{t,x;e) + he 

In this case the energy of the fields depends on e. In fact, we have 

1 



L2{t;e)-Li{t;e) 



gc.,L(y(t,a;)). (52) 



E{t) EE / dx£{t,x) 

Jr ^ 

e{g2{t) ' 9i{t)) 



ay)] +{dyl>{y) 



dy 



L2 — Li 

where we have introduced the "free" fields (the fields when the two mirrors are at rest) 



j=R,L 



= X! / du;ujd^jg^,j{y) + he. 
In the same way, the corresponding quantum Hamiltonian is obtained as 



H{t) = / dxH{t,x) = - 
Jr ' 

e{g2{t) - gi{t)) 



dy 



ay)] +[dyHy) 



L2 — Li 

jdy 

j=l,2 •'^ 



dy dy(j){y) 



L2 — Li 



dy^iy)iy ^ L,) + -<Piy) ] + he 



+ 0(e2). 



(53) 



(54) 



(55) 



(56) 



In this case, the part of the Hamiltonian that describes the interaction between the field and the mirrors is also dependent 
on e. However, in general, it is impossible to adequately describe this part. For that reason, the reactive part of the motion 
force can seldom be calculated. 

For instance, if we consider the generalization to the Lagrangian density (j46|, i.e. 



1 



- {{dt(b)' - id. 



■ Y aj(j)^5{x- Lj{t;e)), 

j=l,2 



(57) 
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then, the part of the quantum Hamiltonian that describes the interaction is 



J = l,2 J = l,2 



e(g2(t)-.gi(t)) 

L2 — Li 



3 = 1-2 



(58) 



And, since 



ay)] +(dyi{y) 



51 - 



is the Hamiltonian of the system when the two mirrors are at rest, that is the "free" Hamiltonian, we can conclude that, in 
the interaction picture, while the mirrors move the full Hamiltonian of the system is given by 



Hi{t) 



e(.92(i)-ffiW) 



L2 — Li 



dy ( dy(t)j{y) ) + X! ( "^/(^j) 



j = l,2 



+ 0(e'). 



(59) 



Finally, we prove that our dissipative part of the motion force coincides with the one obtained in Ill3ll . For times r 
larger than the stopping time, our quantum evolution operator, in the linear approximation, is T"^ = Id. — i J^dt Hi{t). 
We are interested in the term 



(60) 



Integrating by parts, and using the fact that (jjj and are free fields, it follows that 



Ai = -e^ / dtgj{t)Fj{t) + e dt dy 



L2 — Li 



dv<t>i{y) 



2 8^0 7e L2 - Li 



J2 ML,) (dyMh - \5\) - dyMLj + \S\) 



(61) 



Now, from expression (2.4) in Ref. 111211 . 



we obtain that 



lim [dy4>j{Lj + \6\) - dy(l)j{Lj 



^1 = -e ^ / dtgj{t)Fj{t) + e dt dy 

Mt)-9i{t) 



dt- 



L2 — Li 



3 = 1-2 



2a(j)j [Lj 



92{t)-9i{t) 
L2 — Li 
2 



dy(t>i{y) 



(62) 



And, finally, inserting this expression into dtHi{t), we conclude that, for times r beyond the stopping time, we have 

T^=Id + ieY^ 1 dtg-j{t)F-j{t), (63) 

as we wanted to demonstrate. 



3 = 1,2- 
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V. DISCUSSION AND CONCLUSIONS 

By way of the consideration of physically plausible mirrors and of a canonical use of the Hamiltonian approach, we 
have showed that the problem of the negative energies that appears in the Davis-Fulling model can be avoided, the main 
difference with respect to their approach (and other's) being that partially transmitting mirrors, which are transparent to 
sufficiently high frequencies — thus providing a natural and physically sound renormalization — are here being considered. 

We have also proved that our method provides a dissipative radiative-reaction force that fully agrees with the dissipative 
force derived in Refs. ||,27|] and [12]. On the other hand, the motion force calculated using the Hamiltonian approach 
contains some reactive terms, proportional to the the mirrors's acceleration, which do not appear in the results obtained 
in Refs. (2f\ and (l^. Those terms are fundamental in order to ensure that the energy remains positive at any time, 
and consequently, to guarantee the validity of the concept of particle also during the oscillation of the mirrors, which is 
certainly fast but, in the proposed experimental settings, still very small as compared with the speed of Ught lll4ll . 

We have also seen explicitly that albeit a possible (and quite simple) solution to this disagreement could be to perform 
a mass renormalization that completely eliminates the reactive terms proportional to the the mirror's acceleration (see 
Ref. 1IT2I1 ). it would turn out in this case that the definition of particle itself would be impossible to maintain at any time 
while the mirrors move, what would be indeed remarkable, in view of the relatively small velocities involved. Quite on 
the contrary, and in consonance with the realistic boundary conditions imposed by us on the mirrors (which led precisely 
to these additional terms), we take as the most reasonable and physically meaningful renormalization condition to impose, 
to keep those terms in full, in which case the definition of particle during the oscillation of the mirrors can be consistently 
preserved, as well as the fundamental laws of physics too, at any time t during the process. In plain words, as a bonus, the 
fundamental principle of energy conservation holds during the whole evolution towards the end state. 

To finish, we must mention that the dynamical Casimir effect is being discussed right now by many groups in different 
contexts, and that the growing potential of this subject, both as a fundamental phenomenon as well as for the number and 
importance of its applications, is out of question. At very large scales, it is being considered in theoretical Cosmology 
as a most natural explanation of the observed acceleration in the Universe expansion (termed as dark energy) |30]. And 
in a very different context, some laboratory experiments have recently been proposed which would provide an extremely 
nice, alternative proof of the validity of General Relativity and to some semiclassical approaches to Quantum Gravity. In 
addition, they open the way to practical applications of the Casimir effect in nano-electronics and other technologies. In 
those contexts some recent papers have appeared which deserve careful consideration (see e.g. lll4il3l[l ). 

Acknowledgments 

This work has been supported in pait by MEC (Spain), projects MTM2005-07660-C02-01, BFM2006-02842, PR2006- 
0145, FIS2005-25313-E, and by AGAUR (Generalitat de Catalunya), contract 2005SGR-00790. 



[1] V.V. Dodonov and A.B. Klimov, Phys. Rev. A53, 2664 (1996). 

[2] J.Y. Ji, H.H. Jung, J.W. Park, and K.S. Soh, Phys. Rev. A56, 4440 (1997). 

[3] M. Crocce, D.A.R. Dalvit, F.D. Mazzitelli, Phys. Rev. A66, 033811 (2002). 

[4] R. Schutzhold, G. Plunien, and G. Soff, Phys. Rev. A65, 043820 (2002); G. Schaller, R. Schutzhold, G. Plunien, and G. Soff, 

Phys. Rev. A66, 023812 (2002). 
[5] M. Ruser, J. Opt. B7, 100 (2005). 

[6] I. Brevik, K. Milton, S.D. Odintsov and K. Osetrin, Phys. Rev. D62, 064005 (2000). 

[7] S.A. Fulling and PC.W. Davies, Proc. Roy. See. Lond. A348, 393 (1976). 

[8] J. Haro and E. Elizalde, Phys. Rev. Lett. 97, 130401 (2006). 

[9] L.H. Ford and A. Vilenkin, Phys. Rev. D25, 2569 (1982). 
[10] G. Barton and C. Eberlein, Ann. Phys. (NY) 227, 222 (1994); R. Gutig and C. Eberlein, J. Phys. A31, 6819 (1998). 
[11] C. Eberlein, J. Phys. 13, 2151 (1993). 

[12] G. Barton and A. Calogeracos, Ann. Phys. (NY) 238, 227 (1995); A. Calogeracos and G. Barton, Ann. Phys. (NY) 238, 268 
(1995). 

[13] M.-T. Jaekel and S. Reynaud, J. Phys. 12, 149 (1992); 13 1093 (1993); A. Lambrecht, M.-T. Jaekel, and S. Reynaud, Phys. Rev. 
Lett. 77, 615 (1996). 

[14] W.-J. Kim, J.H. Brownell, and R. Onofrio, Phys. Rev. Lett. 96, 200402 (2006). 
[15] M. Razavy and J. Teming, Phys. Rev. D31, 307 (1985). 



14 



[16] C.K. Law, Phys. Rev. A49, 433 (1994); A51, 2537 (1995). 

[17] R. Schutzhold, G. Plunien and G. Soff, Phys. Rev. A57, 2311 (1998). 

[18] H. Johnston and S. Sarkar, J. Phys. A29, 1741 (1996). 

[19] J. Haro, Int. Jour. Theor. Phys. 46, 951 (2007); 46, 1003 (2007). 

[20] M.-T. Jaekel and S. Raynaud, J. Phys. II, 1395 (1991). 

[21] A.A. Grib, S.G. Mamayev, and V.M. Mostepanenko, Vacuum Quantum Effects in Strong Fields (Friedman Laboratory Publishing, 
1994). 

[22] G.T. Moore, J. Math. Phys. 11, 2679 (1970). 

[23] N. Obadia and R. Parentani, Phys. Rev. D64, 044019 (2001). 

[24] N.D. Birrell and C.P.W. Davies, Quantum Fields in Curved Space (Cambridge University Press, 1982). 

[25] E. Elizalde, J. Phys. A34, 3025 (2001); E. Elizalde, M. Bordag, and K. Kirsten, J. Phys. A31, 1743 (1998); M. Bordag, E. 

EUzalde, and K. Kirsten, J. Math. Phys. 37, 895 (1996); M. Bordag, E. Elizalde, K. Kirsten, and S. Leseduarte, Phys. Rev. D56, 

4896 (1997); E. Elizalde, L. Vanzo, and S. Zerbini, Commun. Math. Phys. 194, 613 (1998); E. EUzalde, Commun. Math. Phys. 

198, 83 (1998); Nuovo Cim. 104B, 685 (1989); J. Phys. A22, 931 (1989); E. EUzalde and A. Romeo, Phys. Rev. D40, 436 (1989). 
[26] R.L. Jaffe, Phys. Rev. D72, 021301 (2005); T. Emig, R.L. Jaffe, M. Kardar, and A. Scardicchio, Phys. Rev. Lett. 96, 080403 

(2006); N. Graham, R.L. Jaffe, and H. Weigel, Int. J. Mod. Phys. A17, 846 (2002). 
[27] M.-T. Jaekel and S. Reynaud, Quant. Opt. 4, 39 (1992). 
[28] G.M. Salomone and G. Barton, Phys. Rev. A51, 3506 (1995). 

[29] D.A.R. Dalvit and RA. Maia Neto, Phys. Rev. Lett. 84, 798 (2000); RA. Maia Neto and D.A.R. Dalvit, Phys. Rev. A62, 042103 

(2000). 

[30] E. Elizalde, S. Nojiri, S.D. Odintsov, and P Wang, Phys. Rev. D71, 103504 (2005); E. EUzalde, S. Nojiri, S.D. Odintsov, and S. 

Ogushi, Phys. Rev. D67, 063515 (2003); E. Elizalde, J. Phys. A36, L567 (2003); Phys. Lett. B516, 143 (2001); J. Math. Phys. 35, 

6100 (1994); 35, 3308 (1994); 
[31] F. Capasso, Fellow, J.N. Munday, D. lannuzzi, and H.B. Chan, IEEE J. Sel. Topics Q. Electr. 13, 400 (2007). 



